Automation of Feynman diagram evaluation 
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Abstract 
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On 

^\ ' A C-program DIANA (Diagram ANAlyser) for the automation of Feynman diagram evaluations 

is presented. 

' Recent high precision experiments require, on the side of the theory, high-precision calculations re- 

pL| i suiting in the evaluation of higher loop diagrams in the Standart Model (SM). For specific processes 

thousands of multiloop Feynman diagrams do contribute, and it turns out to be impossible to perform 
these calculations by hand. This makes the request for automation a high-priority task. 

Several different packages have been developed with different areas of applicability. For example, 
FEYNARTS / FEYNCALC @ are MATHEMATICA packages convenient for various aspects of the 
£2 ' calculation of radiative corrections in the SM. There are several FORM packages for evaluating multiloop 

, diagrams, like MINCER and a package Q for the calculation of 3-loop bubble integrals with one non- 

' zero mass. Other packages for automation are GRACE H and COMHPEP §, which partially perform 

full calculations, from the process definition to the cross-section values. 

A somewhat different approach is persued by XLOOPS ||. A graphical user interface makes XLOOPS 
an 'easy-to-handle' program package, but is mainly aimed to the evaluation of single diagrams. To deal 
with thousands of diagrams, it is neccessary to use special techniques like databases and special controlling 
programs. In @ for evaluating more than 11000 diagrams the special database-like program MINOS was 
developed. It calls the relevant FORM programs, waits until they finished, picks up their results and 
QJ , repeats the process without any human interference. 

It seems impossible to develop an universal package, which will be effective for all tasks. It appears 
absolutely necessary that various groups produce their own solutions of handling the problem of automa- 
tion: various ways will be of different efficiency, have different domains of applicability, and last but not 
least, should eventually allow for completely independent checks of the final results. This point of view 
motivated us to seek our own way of automatic evaluation of Feynman diagrams. 

Our first step is dedicated to the automation of the muons two-loop anomalous magnetic moment 
(AMM) \(g — 2) . For this purpose the package TLAMM was developed ||. The algorithm is imple- 
mented as a FORM-based program package. For generating and automatically evaluating any number 
of two-loop self-energy diagrams, a special C-program has been written. This program creates the initial 
FORM-expression for every diagram generated by QGRAF flQ[| , executes the corresponding subroutines 
and sums up the various contributions. In the SM 1832 two-loop diagrams contribute in this case. The 
calculation of the bare diagrams is finished. 

Our aim is to create some universal software tool for piloting the process of generating the source code 
in multi-loop order for analytical or numerical evaluations and to keep the control of the process in general. 
Based on this instrument, we can attempt to build a complete package performing the computation of 
any given process, at least in the framework of a concrete model. 

The project called DIANA (Diagram ANAlyser) jj| for the evaluation of Feynman diagrams is being 
finished by our group at present. 

The program DIANA contains two ingredients: 1. Analyzer of diagrams. 2. Interpreter of a special 
text manipulating language (TM). The TM language is a very simple TeX-like language for creating 
source code and organizing the interactive dialog. 

1 Supported by Bundesministerium fur Forschung und Technologie under PH/05-7BI92P 9. 
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The analyzer reads QGRAF output and passes necessary information to the interpreter. For each 
diagram the interpreter performs the TM-program, producing input for further evaluation of the diagram. 
Thus the program: 

Reads QGRAF output and for each diagram it: 1. Determines the topology, looking for it in the 
table of all known topologies and distributes momenta according to the current topology. If we do not 
yet know all needed topologies, we may use the program to determine missing topologies that occur in 
the process. 2. Creates an internal representation of the diagram in terms of vertices and propagators 
corresponding to the Feynman integrand. 3. Invokes the interpreter to execute the TM-program to insert 
explicitely expressions for the vertices, propagators etc. Producing FORM input e.g., is done by passing 
to it the necessary information from the analyser. 

Executing the TM-program provides apart from the possibility to calculate each diagram using FORM 
or another formulae manipulating language, to do some numerical calculation by means of FORTRAN, 
to create a postscript hie for the picture of the current diagram, etc. 

The program operates as follows: first of all, it reads its configuration hie, which may be produced 
manually or by DIANA as well. This file contains: 1. The information about various settings (file names, 
numbers of external particles, definition of key words, etc.) 2. Momenta distribution for each topology. 
3. Description of the model (i.e., all particles, propagators and vertices). 4. TM-program. Then the 
program starts to read QGRAF output. For each diagram it determines the topology, assigns indices and 
creates the textual representation of the diagram corresponding to the Feynman integrand. All defined 
data (masses of particles, momenta on each lines, etc.) are stored in internal tables, and may be called 
by TM-program operators. At this point DIANA invokes the interpreter and performs the TM-program. 
After that it starts to work with the next diagram. 

When all diagrams are processed, the program may call the interpreter once again (optionally) to 
perform the TM-program a last time. Such last call may be used to do some final operations like 
summing up the results. 

There is a possibility to use the program as an interpreter only. If one specifies in the configuration 

file 

only interpret 

then DIANA will not try to read QGRAF output, but immediately enters the TM - program. 




Figure 1: Topology (-2,2)(-l,l)(5,l)(l,3)(3,4)(4,2)(2,5)(5,6)(3,6)(4,6). 

Topologies are represented in terms of ordered pairs of numbers like (fromvertex,tovertex) (see fig. |l|). 
All external legs have negative numbers. These are supposed to begin with -1. The fist numbers cor- 
respond to ingoing particles, the last numbers to outgoing ones. External legs must be connected with 
vertices of smallest possible identifying number. The number of an internal line corresponds to its position 
in the chain of pairs and the direction from the first to the second number in the pair: 

direction: 5 -> 1 1-^3 3 -> 4 4-^2 2^5 5^6 3^6 4^6 
(-2,2)(-l,l) (5,1) (1,3) (3,4) (4,2) (2,5) (5,6) (3,6) (4,6) 
number: 12345678 



2 



Knowing thus the topology we can assign momenta. Their distribution according to fig. |lj e.g. is added 
like 

topology = (-2, 2) (-1,1) (5,1) (1,3) (3, 4) (4, 2) (2, 5) (5 ,6) (3,6) (4,6) : 
pl,p2,p3,p4,p5,p6,p7,p8; 

This fixes directions and values of all momenta on internal lines. 

The main goal of the TM-language is the creation of text files. In principle we could have used one 
of the existing languages, but we want a very specific language: it should be powerful enough to create 
arbitrary program texts. On the other hand, it should be very simple and easy-in-use, so that even 
non-programmers can use it. 

Similar to the TeX language, all lines without special escape - characters ("\") are simply typed to 
the output file. So, to type "Hello, world!" in the file "hello" we may write down the following program: 

\program 
\setout (hello) 
Hello, world! 

Each word the first character of which is the escape character will be considered as a command. This 
feature makes this language very easy-to-use. 

The user types his FORM program and has the possibility to insert in the same line TM-operators as 
well. For example, the typical part of a TM program looks like follows: 

\setout (d\currentdiagramnumber () . f rm) 

#define dia "\currentdiagramnumber () " 

#define TYPE "\typeO" 

#define COLOR "\colorO" 

#define LINES "\mimberof internal lines () " 

\masses () 

#include def . h 

1 R=\integrand() 

#call f eynmanrules-Q 

#call projectionO 

#call reducing-C'TYPE'} 

#call table-C'TYPE'} 

#call colorfactor{' COLOR'} 

. sort 

drop R; 

g dia'dia' = R; 
. store 

save dia'dia' .sto; 
. end 

\setout (null) 

\system(\ (f orm -1 )d\currentdiagramnumber() .frm) 

Some of the TM - commands are just TM-operators while some are functions (returning a value) 
written in the TM-language itself. For example, \numberof internallines () is a built-in TM operator, 
while the function \masses() is the TM-language function: 

\f unction masses ; \-\let (i , 0) 
\+\do 

#define m\inc(i,l) "\mass (\get (i) ) " 

\while "\numcmp(\get (i) , \numberof internallines () ) " eq "<" loop 
Vend 

This TM-program will generate the FORM input for each diagram and then it will call the FORM 
interpreter by means of the operator \system(). For example, the corresponding part of the FORM 
program generated for diagram number 15 looks like follows: 
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#define dia "15" 
#define TYPE "4" 
#define COLOR "3" 
#define LINES "4" 
#define ml "mmH" 
#define m2 "mmW" 
#define m3 "mmW" 
#define m4 "mmH" 
#include def .h 
1 R= 

l*V(l,mul,mu,2)*(-i_)*em"2/2/s*V(2,0)*(-i_)*l/4*em~2/s~2*mmH/mmW* 
V(3,mu2,+q4-(+q3) , l)=K-i_)*em/2/s*SS(l,0)*i_*VV(2,mul,mu2,+q2,2)*i_* 
SS(3,2)*i_*SS(4,0)*i_; 

#call f eynmanrules-Q 

#call projection-Q 

#call reducing{'TYPE'} 

#call table{'TYPE'} 

#call colorf actor{'C0L0R'} 

. sort 

drop R; 

g dia'dia' = R; 
. store 

save dia'dia' .sto; 
. end 

At present, we have finished the C-part of this project. Also, we have several files with TM-language 
macros to start DIANA. Similar as LaTeX provides the possibility for non-spccialsts to typeset high- 
quality texts using TeX language, these macros permit DIANA to work at very high level. The user can 
specify the model and the process, and DIANA will generate all neseccary files. 
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